Exact thermodynamics of an Extended Hubbard Model of single and paired carriers 
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By exploiting the technique of Sutherland's species, introduced in we derive the exact spec- 
trum and partition function of a ID extended Hubbard model. The model describes a competition 
between dynamics of single carriers and short-radius pairs, as a function of on-site Coulomb repulsion 
{U) and filling (p). We provide the temperature dependence of chemical potential, compressibility, 
local magnetic moment, and specific heat. In particular the latter turns out to exhibit two peaks, 
both related to 'charge' degrees of freedom. Their origin and behavior are analyzed in terms of 
kinetic and potential energy, both across the metal-insulator transition point and in the strong 
coupling regime. 
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I. INTRODUCTION 



In condensed matter, electron systems in regimes of 
high correlation are known to be suitably modeled by 
the Hubbard Hamihonian and its generahzations 
1^-^. For such models, the finite temperature properties 
are the ultimate results which theoretical investigations 
(numerical or analytical) aim to reach, in view of com- 
parisons to experimental data. Indeed some observables 
exhibit intriguing features as a function of the tempera- 
ture, which deserve an accurate interpretation. 
In particular, the thermodynamics of the standard Hub- 
bard model has been widely investigated. In D=l this 
was done by different exact approaches: in pQ,pd[| and in 



for the usual case of nearest-neighbors hopping, while 
|13| for the case of long-range hopping. In dimension 
greater than one recent results were obtained by exact 
diagonalization on small clusters fljjl^ and numerical 
investigations ||l6|,p7t, whereas the case D=c» has been 
examined in |]l8| hy iterated perturbation theory. 
All the results show interesting behaviors as a function 
of temperature, with varying the filling and the Coulomb 
repulsion. This is the case for instance for the specific 
heat, where a double peak structure as well as the ap- 
pearance of quasi universal crossing points were found, 
which features were already noticed in some experimental 
data [[l9y2C| ]. In the strong couphng regime the presence 
of a two-peak structure is usually related to the so called 
'spin' and 'charge' degrees of freedom. Numerical results 
in one and two dimensions ||l^,|l^ show that, at least 
at half-filling, such structure survives also at moderate 
couplings. 

Contrary to the ordinary Hubbard model, which has 
been approached through several techniques, for the ex- 



tended Hubbard models most of finite-temperature re- 
sults have been carried out by means of mean-field theo- 
ries 0. In one dimension, however, it is known that tra- 
ditional approaches to many-body systems such as mean- 
field or Fermi Liquid theories are either unreliable or in- 
applicable. As a consequence, both numerical techniques 
(like Density Matrix Renormalization Group |^l|) and 
non-conventional analytical approaches (like bosoniza- 
tion [p2|) have to be supported by comparison with exact 
solutions, whenever available; this is basically the reason 
for the growing interest devoted to finite temperature ex- 
act results. 

The main technique within exact approaches to one di- 
mensional systems is the Bethe Ansatz (BA), either in 
the Coordinate or in the Algebraic formulation. 
Such technique amounts to guessing for a given model 
eigenstates of the form proposed by Bethe and in 
particular it has been extensively applied to models of 
correlated electrons; for instance, the BA equations for 
a wide class of integrable extended Hubbard models |2^] 
have been recently derived in |^^. However, the actual 
solution of these equations, i.e. the evaluation of the 
quantum numbers characterizing the system (quasi mo- 
menta), is in general quite difficult, and some hypothesis 
on their distribution (string hypothesis |^^) have typi- 
cally to be conjectured. In order to derive the complete 
solution and calculate thermodynamic quantities, one is 
thus reduced to solve a system of infinitely many coupled 
integral equations, which requires dramatic numerical ef- 
fort. More recently, a considerable progress has been 
achieved through the alternative approach of Quantum 
Transfer Matrix , which yields dealing with only a fi- 
nite number of coupled integral equations. This has been 
done for the ordinary Hubbard model ||l^ , for the t — J 
model ISCll, as well as for an extended Hubbard model 
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with bond-charge interaction | |31| . 

Nevertheless, determining the actual properties of a 
model at finite temperature for arbitrary parameters' val- 
ues remains in general a very hard task, even when the 
model is proved to be integrable, and its ground state 
features are possibly derived. 

In the present paper we present the exact thermody- 
namics of a one-dimensional Extended Hubbard Hamilto- 
nian (described in sec. II) whose exact analytical ground 
state properties were obtained in |^ by a technique dif- 
ferent from BA. We called that technique Sutherland 
Species (SS) technique, and here we show how it can 
be exploited to derive explicitly the whole spectrum and 
the partition function of the model (section III). In Sec- 



tion IV we calculate some thermodynamic quantities, 
namely the chemical potential, the compressibility, the 
local magnetic moment, and the specific heat. In partic- 
ular in subsection D we focus on the specific heat, which 
turns out to exhibit a two-peaks structure. The origin 
of such structure and the differences with respect to the 
standard Hubbard model are discussed in section V. 



II. THE MODEL 

The Hamiltonian we are interested in reads: 

(1) 



(1 -n,_^)4^Cj^(l -nj-,,) + Y^cl^cl^Cj 



Here c\^,c^^ are fermionic creation and annihilation op- 
erators on a one-dimensional chain with L sites, tr G {| 
, J,} is the spin label, fija = Cj^Cj^., and (i, j) stands 
for neighboring sites. The Fock space F of the system 
is the product of the L four-dimensional vector spaces 
Vj related to each site j; each Vj is spanned by the ba- 
sis I Dj, I \0)j, I 11) j, which we shall also denote in 
the following as |eQ)j, a = 1, . . . 4 respectively. We shall 
adopt for the ID lattice open boundary conditions; as 
usual, these are not expected to affect the results in the 
thermodynamic limit. 

In the Hamiltonian (|l|) the three terms (which will also 
be denoted as Tit, H-y and 'Hu) represent respectively the 
kinetics of single carriers, the kinetics of paired carriers, 
and the on-site Coulomb repulsion. 

More explicitly. Tit describes the hopping of single elec- 
trons towards empty sites. This term is thus reminiscent 
of the so-called '[/ = cxo Hubbard model'. An important 
difference must be however highlighted: the latter model 
reads VY.{^,J) ,a ^^t^j^V , where V = W.il - h^^fin) 
projects the doubly occupied sites out of the Hilbert 
space (which in that case is actually 3^-dimensional) ; 
in contrast, the term Tit in (0), although not involving 
pairs, does not exclude their presence in the state of the 
system [||. 

The second term in (|l|) is in contrast a kinetic term of 
pairs only; it is worth stressing that the model deals with 
pairs having a very short radius; in fact, while in models 
such as BCS one has several pairs within a radius of the 
coherence length, here the radius of a pair is thought of 
as small with respect to the lattice constant, and is actu- 
ally taken as zero. This kind of term is also used in the 
so called Penson-Kolb-Hubbard model (see |^), where 
one investigates the effects of the pair dynamics without 
explicitly entering the microscopic mechanism yielding 
their formation We also point out that the first and 
the second term in (|l|), though describing the kinetic of 
different kind of carriers (single and pair respectively). 



do not commute at all. 

The third term is traditionally the most important term 
for Hubbard-like models; indeed, according to Hubbard's 
picture, it is the parameter that should drive the metal- 
insulator transition in the d-transition metal compounds. 
Loosely speaking, the ratio U /t can be thought as pro- 
portional to the inverse of the pressure applied on the 
sample: by increasing the pressure one reduces the lat- 
tice spacing and thus makes the hopping amplitude more 
relevant with respect to U . 

The first two terms of the Hamiltonian are in general 
competing: indeed Tit would favor delocalized waves of 
single carriers, avoiding the formation of pairs; Tiy low- 
ers instead the energy when electrons form tightly bound 
pairs moving along the chain. This competition is in ad- 
dition modulated by both the term in U and the filling, 
i.e. the density p of electrons in the chain. This can be 
seen by examining the case 

Y = -t . (2) 

Indeed for this value of the coupling constant the model 
has been proved to be integrable and the exact 
ground state phase diagram (reported in fig. ^ has been 
obtained in 0|. Tuning U and p the model exhibits in- 
teresting features; for instance, even when the value of 
filling is p < 1 and at moderate {U < 2t) Coulomb repul- 
sion, it is energetically favorable for the system to form 
pairs and let them move instead of having only singly 
occupied sites. 

In region I the ground state (GS) is made of only doubly 
occupied and empty sites; in region II we have also singly 
occupied sites (either | "f) or | J,)). In region Ill-a the GS 
is that of the U = oo Hubbard model and is made of 
singly occupied sites (metal). In region Ill-b the GS of 
the model reduces to that of the atomic limit of the Hub- 
bard model (insulator). At half-filling {p ~ 1) a charge 
gap Ac = [/ — 2t opens for any U > 2t. 
We wish to stress that, unlike many exactly solved elec- 
tron systems, the model (IT]) is not particle-hole invariant: 
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indeed the first term breaks up the invariance; this leads 
to the shape of the phase diagram shown in fig.0, which 
is asymmetrical with respect to half- filling. 




P (filling) 

FIG. 1. Ground state phase diagram of the model for 
Y = t, from [1]. Empty, bared, and full circles respectively 
represent empty, singly occupied, and doubly occupied sites 
in the ground state. 



III. SPECTRUM OF THE SYSTEM 

In the following we shall assume Y = —t, since such 
relation allows for the integrability, as observed above. 
In this case, the Hamiltonian (|l|) can be rewritten in the 
form 



(3) 



(ij) 



where Tij accounts for the first two interaction terms 
in (|l|). The term in U is easily checked to commute 
with j) "^i-i- Due to the condition (|^) Ti.j, exhibits 
the structure of a generalized off-diagonal pcrmutator be- 
tween Physical Species (PS), which arc the local vectors 
|eQ)'s. More explicitly, while the ordinary off-diagonal 
permutator, when acting onto \ea)j ® |e^)j-(-i, returns 
\e[3)j ® \ea)j+i for any a ^ (3 and zero for a = P, a gen- 
eralized one makes the exchange or gives zero according 
to the specific values of a and p. In our case, Tij per- 
mutes the PS of two neighboring sites only if one belongs 
to the group A and the other to the group B, where 



A = |T),U),UT) , B = \o) 



(4) 



In all the remaining cases T^j gives zero. The above 
groups A and B of PS can be identified with the Suther- 
land Species (SS) of the model (|l|)(see 0); the notion of 
SS is strictly related to the structure of the Hamiltonian 
and not to that of the underlying Hilbert space . 
In D=l a generalized permutator between PS has the 
same eigenvalues as an ordinary permutator between the 
corresponding SS. This is actually what allows to provide 
the exact spectrum, as we shall see below. 



The Fock space F of the system is F = ®^=ol^^' 
where is the A^-electron Hilbert space {N = 

J2i=i "-it + "-ii)- However, due to the properties of the 
Hamiltonian, it turns out to be useful to rearrange F in 
terms of Hjv^, i.e. the spaces spanned by all the vec- 
tors that have a definite number Na of sites occupied 
by states of species A ('A-sites' henceforth). Clearly 
Nb = L — Na- According to the properties of gener- 
alized permutator fulfilled by Ti, the latter commutes 

with Na = SiLi "-it + "ii ~ f^i^f^ii-i Etnd thus H^v^ is 
preserved by the dynamics (this would hold in any di- 
mension). In addition, dealing with an open chain, one 
can have 3^^ possible sequences S of A-sites for a fixed 
number Na- Notice also that, since i) the first term of (j^) 
only permutes A with B and gives zero otherwise, and 
ii) the second term merely counts the number of doubly 
occupied sites, also the sequence S is preserved by the 
dynamics, and it can be identified with an invariant sub- 
space within H Ma ■ The dimension of each of these 3^-^ 
subspaces is (j^^), accounting for all the possible actual 
positions of the A-sites along the chain. One can repeat 
the above foliation for all the H^Va's {Na runs from to 
L) and check that the Fock space is completely recovered: 



L 

Na 



= 4^ 



(5) 



so that F = 0^^^o 

Focusing on a given H i ons can characterize each of its 
basis vector by specifying two discrete-valued functions: 
S{m) and J(m) (m = 1, . . . Na)- The former, which is 
valued 1 (for | t)), 2 (for | [)) or 3 (for | If)) determines 
the sequence S of A-sites, and thus the invariant sub- 
space in which the vector lies; the latter, which is valued 
1 to L, determines the actual positions of the m-th A- 
site along the chain. The basis vectors can therefore be 
referred to as \{S}, {J} ), where "{ }" is to remind that 
S and J are functions. 

In realizing that the Hamiltonian can be separately diag- 
onalized within each subspace characterized by a given 
v4-sequence 5, it is also crucial to observe that each of 
such invariant subspaces can be put in a one-to-one cor- 
respondence with the states of Na spinless fermion space 
(or equivalently with a spin 1/2 model with magnetiza- 
tion L — Na) as follows 



{S},{J}) 



' Na 

n 



'Jim) 



|0) 



(6) 



where are the creation operators for a spinless fermions 
and {S} the sequence of the subspace. 
Similarly to what has been done in Q for another ex- 
tended Hubbard model, it is also easy to derive the form 
of an effective Hamiltonian for the spinless fermion states: 
indeed, since the first term in (^ reduces to a permutator 
between SS, it actually acts on the considered subspace 
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in the same way as a free Hamiltonian — i j) ^l^j ^'^^^ 
on the spinless problem space. The second term simply 
counts the number of species A of kind | jt); namely 
N-^l = X^i^i "'iT'^ii = N — Na- Therefore the spectrum 
in each subspace is given by: 



i;({n^};A^)= ^(-2t cos fc - ?/) ; 



UN 



(7) 



fc=i 



where {n^} are quantum numbers valued or 1, fc = 
ttI / {L + 1) {I = 1, . . . L), and N is the total number of 
electrons (which ranges from Na to 2Na)- The eigen- 
vectors are given by the anti-transformed through (||) 
of spinless fermion eigenstates (Hfc Si^Li sin(fcz)a|)|0), 
where the product is over Na of the L allowed values of k. 
When passing from a subspace of \~\na to another, one 
finds identical replica of this spectrum, which amount 
to having a degeneracy of the eigenvalues. The degen- 
eracy g corresponds to the different ways in which one 
can choose a species A at a given site provided that N 
remains unchanged (i.e. one has the freedom to change 
singly occupied | |) into | I) and viceversa); it is therefore 
easily seen that 



2Na 



{"fc } ^ = ^A 

L 



2NA-N I Na 

N - Na 



^ (2 + e-'3(c^-A'))JV.4 ^ [I3{2t cos k + 



n,. 



2Na-N 



Na 
N ~Na 



(8) 



To conclude this section, we wish to emphasize that the 
spectrum (^) has been derived by means of Sutherland 
species technique under open boundary conditions. In 
fact the same model was also studied under periodic 
boundary conditions |2^], within the Algebraic Bethe 
Ansatz approach. However in the latter case the result- 
ing equations for the quantum numbers do not allow a 
straightforward evaluation of the eigenvalues; indeed the 
thermodynamics of had not been derived yet. 



IV. THERMODYNAMICS 



Thanks to the exact spectrum obtained in the previ- 
ous section, we can now pass to the study of its ther- 
modynamics, through the exact calculation of the grand 
partition function. The language of Sutherland's species 
turns out to be very useful to this aim; indeed, due to 
the rearrangement of the Fock space described above, one 
can write 



exp —(3 



^(-2icos/c - U): 



.k=l 



fi{U~ij)N 



W [l + ex:p[(3 {2tcosk + II + v{U, 13, pi))]) 



(9) 



fe=i 



where we have defined v{U, f3, p) = ln(2 + e^''^^"^))//?, 
(3 = l/{kBT) being the inverse temperature and /i the 
chemical potential as usual. 



The grand-potential is easily obtained as uj ^ 
Lo{(3; U; fi) = — limi^oo (3~^{lnZ / L). After introducing 
pLeff ~ p + V, OJ reads 



1 r 

uj{P;U;p) = -/ dk \n{l+ eiip[p{2tcosk + peffiU,(3,p))]) (10) 



Remarkably, the grand potential is formally similar to 
that of a tight-binding model with an effective chemi- 
cal potential Peff- We stress that p^f f{U, [3, p) depends 
on the on-site Coulomb repulsion, the temperature and 
the chemical potential, in a highly non-linear way. This 
yields peculiar features of the model, as we shall show in 
the following. 

In deriving the thermodynamics of the system, it is cus- 
tomary to eliminate p in favor of the filling p; the latter 
can be computed as p = —duj/dp, and the result turns 
out to be of the following form 

p(C/,/3,m) = {l + C{U,P,p))-pA{(3,Peff{U,l3,p)) (11) 
where: 

exp (-/?([/ -A^)) 



C{U,(3,p) 



2 + eyi^{-(3{U ~ p)) 



(12) 



and 



PA{(3,Peff) 



dk 



1 -I- exp [/3(— 2t cosfc — Peff)] 



(13) 



Notice that differentiating uj with respect to pef / instead 
of to p would yield only the right factor pA of (|lT|); the 
non-linearity of as a function of p results in the ap- 
pearance of C in the left factor; this causes the relation 
p = p{p; T; U) implicitly defined by ( [ll| ) to be very dif- 
ferent from that of a tight binding model, as we shall 
explicitly show in next section. 

The two factors in ( pi] ) deserve some comment; pA is 
nothing but the density of A-sites along the chain, defined 
as PA = ^vaiL^+oo{N a) I L] the functional dependence of 
PA on (3 and /Xe/ / is that of a spinless tight-binding model. 
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The left factor provides information, through the func- 
tion C, about the kind of occupancy of the sites of the 
chain; indeed when C ^ most of the occupied sites are 
singly (s) occupied, whereas if C ^ 1 most of the occu- 
pied sites are doubly (d) occupied; intermediate values 



Eq.(^^ naturally allows to identify in £ a kinetic en- 
ergy IC and a potential energy V. The former is defined 
as the weighted integral of — 2icosfc and the latter as 
the weighted integral of —U, which actually gives ~UpA, 
according to (^3|). In fact the actual potential energy 
would also contain the last term Up of (^4|); however, 
since this is merely a constant with respect to tempera- 
ture, we prefer not to include it in the definition of V, so 
that the latter describes the only temperature dependent 
part of the potential term Uhi^fiii. Notice that with this 
choice the potential energy is attractive for positive U . 
Notice also that, although /C and V are clearly related to 
the hopping terms and to the on-site Coulomb repulsion 
respectively, they are not mutually independent: indeed 
/C depends not only on t but also U and viceversa for 
v. We shall come back to this point in discussing the 
specific heat in section V. 

A. Chemical potential 

The chemical potential /i(p; T; U) of our model is 
shown in fig.^ aX U ^ t (a) and U ~ At (b), for dif- 
ferent values of the temperature. 

Focussing first on the solid curves, representing the case 
T = 0, one can realize that even in the ground state 
the relation between /i and p is quite different from 
that of a spinless tight-binding model, which would read 
p-{p;T = 0) = -2tcos(7rp). 

In particular, in fig. (a) we notice that a 'plateaux' ap- 
pears, in correspondence with region II of the ground 
state phase diagram (see fig.0) . Interestingly, such shape 
reminds that of a coexistence region, connecting the 
phase of single carriers (region Ill-a) to that of pair- 
carriers (region I); this would imply that, as the filling 
is increased, the ID lattice starts exhibiting macroscopic 
regions made of only single carriers separated by other 
macroscopic regions where only pairs are present. In 
fact, eigenstates with such features are certainly present; 
however, they are degenerate with other eigenstates, in 
which single and pair carriers alternate with no macro- 
scopic order. This is basically due to the degeneracy of 
^-sequences in such region. 

In fig.(b) a vertical jump is instead present at half-filling, 
as an hallmark of the opening of the charge-gap. The flat 
part of the solid curve for p > 1 just coincides with the 
atomic limit behavior (region Ill-b of fig.|l|). 
Considering now the curves at finite temperature of figJ|, 



indicate the percentage of d with respect to s sites. 
To conclude this section we wish to comment about the 
energy (per site) of the system; the latter is obtained by 
£ = — limL^oo d{lnZ /L)/d/3 + pp and reads 



(14) 



one can observe how the edges present at T = smoothen 
as soon as T > 0. A remarkable feature is the presence 
in fig.(b) of a nearly universal point {p* = 4/3, /i* = U), 
where all the curves of sufficiently low temperatures ba- 
sically intersect. Such kind of points are in general 
determined through the conditions dp/dT — and 
d^p/dT^ ~ 0. It is in fact possible to calculate that 
for any U > 2t and p > I (region Ill-b of fig.^ the 
low-temperature behavior of p is given by 

p^U + kBT\ni-^^—-l\+0{e- — ) , (15) 

whence the above conditions are both fulfilled up to ex- 
ponentially small terms in ksT/t. 

We shall also see in subsection D that nearly-universal 
crossing points are exhibited by other observables of the 
model, such as the specific heat. 




/ '; 

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



P (filling) p (filiing) 

FIG. 2. The relation between p and ^ for different temper- 
atures, at [//t = 1 (a), and U/t = 4 (b). For T = 0, the curve 
in fig. (a) shows a plateaux, related to 'phase coexistence' in 
mixed region II of fig. whereas the curve in fig.(b) exhibits 
a jump in p, due to the opening of charge gap at half-filling. 

Eq. ( p^ also points out that in our model a linear low- 
temperature behavior is possible, differently from the 
tight-binding model, where only even powers in T are 
allowed in the Sommerfeld expansion. In general, in our 
model, different behaviors of p arise according to the val- 
ues of U and p. For instance, for U and p belonging to 
the mixed region II of figj^, the chemical potential at 
low-temperature has again a linear term. 



£(U,(3,p) = - 

TT 



dk 



'2t cos k ~ U 



1 -1- exp I 



-2tcosk- iieff{U,(3,p))] 



Up 
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M -CZ + fesT In (^^|-^^ +0((A;BrA)') , (16) 

but with a coefficient which depends on C/, since p = 
p{U) ^ Ti-^ cos-^{-U /2t). 

In contrast, when the charge gap Ac = C/ — 2i opens {i.e. 
at p = 1 and U > 2t), p, acquires a highly non-hnear form 



(17) 



indicating that the behavior is definitely different to that 
of an intrinsic semiconductor. 



singularity is indeed present for all values of U and p 
that belong to that region of the ground state. The di- 
vergence of K can be proved to be of the type oc T^^. 
In contrast the behavior for T ^ at p = 0.5 and p = 1.5 
is regular since such filling values belong to regions Ill-a 
and I respectively. 




1.2 1.6 



k^T/t 

FIG. 3. The behavior of /i as a function of temperature, for 
fixed U/t — 1 and different values of filling. The asymmetry 
with respect to the half-filled case p = 1 is ascribed to the 
lack of particle-hole invariance of the model. 

In fig.|^ we explicitly examine the behavior of p as a func- 
tion of temperature, for a fixed value of on-site Coulomb 
repulsion {U/t ~ 1) and for different fillings. A main dif- 
ference has to be emphasized with respect to the case of 
a tight-binding model: in the latter the curves of p are 
specular for filling values that are symmetric with respect 
to half filling (i.e. p —p for p — > 2 — p), whereas this 
is not the case in our model, due to the fact that it is not 
particle- hole invariant. 




p = OS 



1.0 

k„T/t 



FIG. 4. The compressibility as a function of temperature, 
at fixed on-site Coulomb repulsion U/t = 1.0 and for different 
filling values, k diverges as T — > for values of U /t and p 
belonging to mixed region II of fig. |^. 



B. Compressibility 

The compressibility k — dp/ dp can be easily evaluated 
through eq.(|ll]). In fig.^ we have plotted func- 
tion of the temperature, for a fixed value of U (namely 
U/t ~ 1.0) and for different fillings. One can observe the 
change in the low-temperature behavior when tuning the 
filling: at p = 0.5 the behavior is regular, while at half- 
filling K undergoes a singularity for T — > 0; eventually 
(p = 1.5) its behavior is again regular. The reason for 
the low-temperature singularity at p = 1 is that in the 
ground state the point {U/t = l:p — 1) is situated in 
region II (see fig.^, i.e. in the region where the chemical 
potential exhibits the plateaux, as shown in fig.H such 



In fig.^ we have examined in detail the case of half- 
filling, plotting function of T for different values 
of U; one can explicitly observe how U = 2t is the crit- 
ical value separating the divergent behavior for U < 2t 
from the regular one for U > 2t. Indeed, as soon as 
U > 2t, the divergence becomes a pronounced peak in k; 
the temperature T* at which the peak occurs increases 
with increasing J7, similarly to what happens in the ordi- 
nary Hubbard model, according to the results jlj]. No- 
tice that in contrast no singular behavior is expected at 
moderate U' s in the ordinary Hubbard model at half- 
filling, since in that case the system is insulating for any 
positive U . 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



k^T/t 

FIG. 5. The compressibility as a function of temperature, 
at fixed filling p = 1 and for several values of on-site Coulomb 
repulsion. As the charge-gap opens {U > 2t) k acquires an 
exponential low-temperature behavior. 



C. Local magnetic moment 



respect to the temperature. Passing through U = 2t, 
an abrupt change in the low-temperature slope occurs: 
the curves of Ao suddenly flattens as soon as U > 2t. 
This reflects the metal-insulator transition occurring in 
the ground state; indeed the opening of the charge gap 
causes the formation of pairs to be highly unfavored at 
low T's. 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

k„T/t 



The local magnetic moment was first introduced in fl^ ] 
and is defined as 



Ao = lim {j - TT-jlf) ■ 



(18) 



FIG. 6. Local magnetic moment as a function of the tem- 
perature, at half filling for several values of U /t. Notice how 
the low-temperature behavior abruptly changes across the 
metal-insulator transition point. The figure indirectly pro- 
vides also the behavior of pA, since Ao = 2pA — p. 



It characterizes the magnitude of spin at each site, i.e. 
the degree of localization of electrons. In terms of density 
of A-site, Ao can be easily rewritten as Aq = p — 2p|| — 
2pA — where pA can be computed from eg. (|l^) . In fig.^ 
we have reported the local magnetic moment at half fill- 
ing for different values of the on-site Coulomb repulsion. 
One can observe that the behavior of Aq, even within 
a relatively small range of values of C/, is quite rich. In 
order to describe it, we first consider the case of small val- 
ues of U (namely U = l.At in the figure); we recall that 
in the ground state such value corresponds to the mixed 
region II (see fig. I at p = 1), meaning that hopping 
paired electrons are present at T = 0; as the temper- 
ature is turned on, Aq first increases with T (indicating 
that the pairs are broken in favor of single carriers); how- 
ever, after reaching a maximum at a temperature T*, Ao 
starts decreasing for higher T's, denoting that pairs are 
now reformed by higher thermal excitations. Accord- 
ing to the above observations, it easy to realize that the 
temperature T* decreases with increasing U ; in fact the 
maximum disappears for U ~ 1.85, so that Aq becomes 
a definitely decreasing function of the temperature. At 
U = 2t, Aq reaches at T = the saturation value 1 (all 
singly occupied sites), with an infinite derivative with 



D. Specific heat 

In this subsection we present our results on the specific 
heat of model (pi) which can be computed through 



_ d£ _ ^ ( d£ d£ dp dp 



(19) 



where the energy £ is given by (|lj). Below we study the 
temperature dependence of Cy when varying the physical 
parameters U and p. The exact calculation shows that in 
our model a two-peak structure is definitely present not 
only in the strong coupling regime, but also at moderate 
C/'s. 

We start considering the case of half-filling (p = 1). The 
two peaks appear first for 1.3 < U/t < 1.8 (see fig.0); 
in this range of ?7's, Cy also exhibits a nearly univer- 
sal crossing point at ksT 0.85t; we shall comment on 
such feature at the end of this subsection. The peaks 
eventually merge into one for U/t ^ 1.85. However, as 
soon as U > 2t (see fig.^), a new well pronounced low- 
temperature peak appears. The recovered double-peak 
structure is present up to J7 ~ 3i, where finally only one 
peak survives. 
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By comparing figs.|7| and|^, one can notice that the metal- 
insulator transition point U = 2tis also the hallmark of a 
crossover in the low-temperature behavior of Cy ■ In par- 
ticular, the calculation shows that for U < 2t the latter 
is linear, 



(20) 



where p is defined as in eg. (Jig) . In contrast, for U > 2t, 
Cv exhibits an exponential-like behavior given by 



where Ac = [/ — 2Hs the charge gap. 
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FIG. 7. Specific heat as a function of T at half filling, for 
different values oi U/t below the metal-insulator transition 
value: a two-peak structure is present, as well as a nearly 
universal crossing point. The low-temperature behavior is lin- 
ear. 
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FIG. 8. Specific heat as a function of T at half filling and 
different values of U /t, just above the metal-insulator transi- 
tion point: the two-peaks structure definitely disappears for 
U > 3t. The low-temperature behavior is exponential for 
U > 2t. 



To conclude the study at half-filling we have exam- 
ined the case of large U/t (see figj|). The result shows 
that only one peak is present, at a temperature which 
increases almost linearly with U {ksT ^ 0.21C/). This 
result can be understood considering that at large U/t 
the spectrum (Q) exhibits two different energy scales: i) 
a low-energy scale (~ t), which describes fluctuations in 
the A-band, whose effective filling is given by the value 
of PA] ii) a high-energy scale (of the order of U) involv- 
ing the formation of on-site pairs, favoring the decrease 
of the number of A-sites. The former channel is actu- 
ally active only for p < 1, since at half- filling the A-band 
becomes completely filled: indeed in this case we have 
PA — i for fcsT ~ t, as can be deduced from fig. ^ of the 
local magnetic moment at large U/t. 
Only the high-energy channel is thus active, and its con- 
tribution is well described by the atomic-limit model (i.e. 
t — Y — 0), shown by the dotted curve in fig.^ The 
slight deviations are due to the fact that, as pairs are 
formed from singly occupied sites via thermal fluctua- 
tions, the number of effective species A decreases, and 
the formed A-holes can produce (relatively small) fluctu- 
ations with T. However, the larger is U/t, the better is 
the agreement with the specific heat of the atomic limit. 
We also wish to emphasize that the behavior is differ- 
ent from that of the ordinary Hubbard model, where two 
peaks appears at low temperatures in the strong coupling 
limit at half filling. In fact, although in the Hubbard 
model the lower Hubbard band is filled, spin excitations 
of low-energy (~ J = At'^/U) are active. These kind 
of excitations are instead absent in our model; we shall 
comment in more detail in section V about this point. 
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FIG. 9. Specific heat as a function of T for fialf-filling and 
4:t < U < 16t. In the strong coupling regime, the two peaks 
structure disappears: the remaining peak is well described 
from the atomic limit model (dotted curve). A similar behav- 
ior is obtained also at any p > 1. 

In fig.^ we investigate the specific heat for filling val- 
ues below half-filling, namely p = .75. 
As fig.|l^(a) shows, a double peak structure of Cy ap- 
pears; however, two important differences have to be em- 
phasized with respect to the case of half-filling: in the 
first instance, here the double-peak structure arises and 
becomes more evident for large values of C/'s, whereas at 
half filling it is present at moderate [/'s; secondly, the 
temperatures of the two peaks are quite higher than the 
corresponding ones of the half filled case. In particu- 
lar the position of the low-temperature peak is practi- 
cally independent of U, whereas the high-temperature 
one strongly depends on it, similarly to what happens 
for the only peak present at half-filling in the strong cou- 
pling regime (fig]^)- 

The two peaks of fig. (a) have to be related to the two 
energy scales emerging in the spectrum when U ^ t, as 
discussed above; in particular, the low-temperature one 
is attributed to the fluctuations of the A-band, which is 
now partially filled, unlike for half-filling. We recall that 
in this range of the parameters U and p, the ground state 
of the model is that of the U — oo model (region Ill-a 
of fig.^; since the formation of pairs is strongly inhib- 
ited for high [/'s, the physics of low-energy excitations 
is fairly captured by that of the U — oo model at finite 
temperature, as shown by the solid curve in fig. (a). In 
flg.(b) the case U/t — 8 is examined in detail; in this case 
the sum of the specific heats of [/ = oo model and atomic 
model practically recovers the actual Cy of our model. 
Such agreement improves with increasing [/, whereas at 
moderate values of U the argument of energy scale sepa- 
ration does not hold: indeed the high temperature peak 
merge into the low-temperature one for U ^ 2t, and Cv 
is no more given as the sum oiU = oo and atomic limits 
(see fig. (c)). 
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FIG. 10. The specific heat as a function of temperature for 
p = 0.75. In fig. (a) Cv is plotted for diflterent values of U/t; 
(b) at strong coupling {U /t — 8) the specific heat of the model 
(solid line) is well reproduced by the sum (dot-dashed) of the 
specific heat of U = oo model (dashed) and of the atomic 
limit (dotted); (c) this is not the case at moderate coupling, 
where the energy scales of the two models become comparable 
{U/t ^2). 

We have also considered the case of filling values 
greater than one. In the strong coupling regime the 
ground state has the A-band completely filled, the sites 
of the chain being all occupied (either singly or doubly, 
as shown in region Ill-b of fig]|); the low -energy scale is 
thus frozen, just like in the case of half-filling. This yields 
the specific heat behavior be actually described by that 
of the atomic limit, similarly to fig.|9|. The temperature 
of the peak grows linearly with U {ksT ~ c{p)U), the 
coefficient c being an increasing function of the filling p. 

Figs.|ll| and |l^ examine the filling dependence of the 
specific heat at fixed coupling values. More precisely, fig. 
|TT| reports the results obtained in the strong coupling 
case. As anticipated above, in this case the low tempera- 
ture peak is perfectly recovered from the U — oo model; 
notice that, since the latter is particle- hole symmetric 
around quarter filling (p = 0.5) the low-temperature be- 
havior of curves related to filling values that are sym- 
metric with respect to p = 0.5 is basically identical. In 
contrast, the higher temperature peak does not exhibit 
such symmetry, being related to the atomic limit of the 
Hubbard model, which is no more particle-hole symmet- 
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ric around quarter filling. 

Fig.|l^ is concerned with the behavior at moderate J7's 
(namely U/t — 1.5) as a function of p; the remarkable 
feature is the appearance of a nearly universal crossing 
point at low temperature (ksT 0.2 1), for a finite range 
of filling values (1.0 < p < 1.3). Similarly, a nearly uni- 
versal crossingpoint also occurs at fixed filling for vary- 
ing U, as fi-g|3 shows. The latter type of behavior is 
also e xhibite d by the ordinary half-filled Hubbard model 

T2i,r 



3 17 



16 



however, to the authors' knowledge, theo- 
retical investigations were mostly limited to the case of 
fixed filling and varying U/t. In contrast, here we have 
explored the case of varying p as well; this is interesting in 
view of comparison with experimental results, where U/t 
can be roughly interpreted as the inverse pressure, and 
6 = \1 — p\ as the doping. In fact, this type of universal 
behavior has been observed in many heavy-fermion com- 
pounds, such as cerium ones, both at fixed doping with 
varying pressure , and at fixed pressure with varying 
doping j2^. Let us notice that, for the ordinary Hub- 
bard model, the presence of the nearly universal point 
in U has been explained in , as a consequence of the 
fact that the entropy S at high temperatures does not 
depend on U, in that case. For our model, S at high 
temperatures is also independent of [/; however it turns 
out that it does depend on p. Hence we expect that the 
argument in cannot be applied to explain the nearly 
universal crossing point in p shown in fig.[l^. 



p = 0.9 




FIG. 11. The specific heat as a function of temperature at 
strong coupling and p < 1. The low-temperature behavior is 
the same for values of p symmetric with respect to p = 0.5: 
indeed in this case low-energy excitations are well described 
by the U = oo Hubbard model, which is particle-hole invari- 
ant around quarter filling. Differences instead emerge at high 
temperatures. 



p = 0.8 




FIG. 12. The specific heat as a function of temperature in 
the moderate C/'s regime. A nearly universal crossing point 
with varying p at fixed U is observed for values of p in the 
range 0.9 < p < 1.3. 




FIG. 13. The specific heat as a function of temperature at 
half filling for negative values of U. Fig. (a): the double peak 
emerges as |?7| is increased. Fig.(b): in the strong coupling 
regime {U/t = —15.0) Cv is fairly reproduced by the sum 
of the XXO contribution (dashed line) and the atomic limit 
contribution (dotted line). 

Finally, the specific heat Cy is investigated in fig.p^ 
also for negative values of the Coulomb interaction, at 
half-filling. The behavior is quite different with respect 
to the positive U case for moderate and intermediate U 
values, since no double peak is present. 
In contrast, such structure emerges at higher coupling 
values; also in this case two separate energy scales 
emerge. However, the low temperature peak is now re- 
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produced by that of the XXQ model {t = 0), whose 
ground state actually coincides with that of our model, 
for these values of U and p (see region I in fig.|l|). The 
high-temperature peak is still due to the ncgative-C/ 
atomic limit {t — Y = 0). In figure (b) it is clearly 
shown how, in the strong coupling case, the simple sum 
of the specific heats of XXO and atomic limit perfectly 
reproduces the result for our model; this is not the case 
by at lower U values. 



V. DISCUSSION 



As outlined in the previous section, our results show 
that the specific heat exhibits a two-peak structure for 
different values of on-site Coulomb repulsion U and fill- 
ing p. In the present section we wish to discuss the origin 
of the two peaks, since in last years much effort has been 
made to clarify a similar behavior occurring in the ordi- 
nary Hubbard model. As mentioned in the introduction, 
in the latter model the two peaks are usually explained 
in terms of 'spin' and 'charge' excitations. 
The above argument cannot be applied here, since our 
model involves only 'charge' degrees of freedom: in fact, 
from the formal point of view of quantum numbers n^, 
the excitation processes in the spectrum have the 
typical feature of 'charge' excitations (in the sense of 
A-species). It is however worth emphasizing that, just 
like for the ordinary Hubbard model, the nomenclature 
based on quantum numbers does not strictly correspond 
to its physical meaning. In our case, the 'charge' degrees 
of freedom of A-species actually carry both charge and 
spin density fluctuations: the break-up of a localized pair 
into two single carriers indeed leads to a redistribution of 
the charge density as well as to the formation of a triplet 
replacing a singlet state. 

In our model any peak of the specific heat has thus to 
be ascribed just to 'charge' excitations. We have seen in 
section IV-D that, when varying the parameters U and 
p, the peaks can merge into one, and possibly reappear. 
In the following we shall discuss such a structured be- 
havior through the kinetic and potential contributions to 
Cv, namely the derivatives /C' and V' with respect to the 
temperature of IC and V, defined when giving the inter- 
nal energy (|lj). 

We start by the case of the strong coupling {U 3> t), 
where our results show a two-peak structure for positive 
U and p < I (see figTO-(b)), as well as for negative U 
at any filling (see fig.|l3|-(b)). Since in these regimes the 
characteristic energy scales of the kinetic term (t) and the 
potential term (J7) of the Hamiltonian are well separated, 
it is expected that each of the two peaks is related to one 
of these terms. In fig.|l^ we have thus plotted /C' and V' 
for U/t = 16 and p = 0.75: the two peaks are indeed in 
perfect correspondence with the contributions of /C and 
V. It is also worth stressing that these two contributions 



can be quite well described at strong coupling in terms 
of two different models: explicitly, the low-temperature 
kinetic behavior is captured by the U = oo model for pos- 
itive £/'s (fig. 
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-(b)) and by the XXO model for negative 
U's (fig.|13|-(b)); the high-temperature potential behavior 
is instead described by the atomic limit. 




k,T/t 

FIG. 14. The kinetic (dashed) and potential (dotted) con- 
tributions to the specific heat (solid), at strong coupling 
{U/t = 16), for p = 0.75. The low-temperature peak is basi- 
cally due to IC' , while the high-temperature peak stems from 
V' ■ In this regime {U ^ t), K! is also well described by the 
specific heat of 17 = oo Hubbard model, and V' by that of the 
atomic limit (dashed and dotted curve of fig.|lo|). 

In contrast, in the regime of moderate C/'s the two en- 
ergy scales become comparable, and the above argument 
is not applicable. This gives rise to a completely differ- 
ent scenario; for instance, at half-filling we observe that 
by lowering V the single strong-coupling peak splits into 
two, whereas for p < 1 the two strong-coupling peaks 
merge into a single one. In practice, while for |?7| ^ t 
the kinetic and potential terms decouple, at moderate 
t/'s it is the competition between the two kinds of energy 
that determines the actual shape of the specific heat. 
This can be understood by recalling the structure of 
the energy spectrum (see eq.(|^)); both terms can be ex- 
pressed in terms of the quantum numbers n^, where the 
total number of A-sites is not a fixed quantity, but can 
vary in the range Na € [A^/2;Af] (the electron number 
N being obviously fixed). This property actually yields 
the competition between V and /C: indeed the kinetic 
term may favor the decrease of X A^ in order to eliminate 
possible positive contributions of — 2tcos/c, whereas the 
potential term favors the increase of X a (*-e. the break- 
ing of on-site pairs). This competition is already active 
at r = 0, causing the appearance of the different regions 
in the ground state phase diagram. 

At finite temperature two more mechanisms enter driving 
such competition: i) the density pA of A-carriers varies 
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with r, according also to the values of U and p; ii) the 
kinetic term exhibits the usual thermal fluctuations. The 
former represent the crucial difference with respect to an 
ordinary free spinless fermion model, where only thermal 
excitations are present, at fixed number of carriers. No- 
tice also that the variability of pA can happen to contrast 
the effect of thermal fluctuations: this is the case when 
Pa decreases with T, since this would yield a reduction of 
/C, while thermal fluctuations would lead to an increase 
of it. As a consequence, a further competition, concerned 
with the purely kinetic contribution, may occur. 
In fig.|l5| we plot the derivatives K,' and V' of the kinetic 
and the potential parts for various moderate U^s at half- 
filling. Starting from U/t = 1.6 we observe that at low 
temperatures both IC' and V' exhibit a peak at nearly 
the same temperature Ti; this is due to the fact that 
in this regime they are driven by the same mechanism 
(formation of pairs from singly occupied sites). The two 
contributions of opposite signs do not completely can- 
cel each other; the kinetic one prevailing, a kinetic low- 
temperature peak appears in Cy- Notice that the value 
of Cv at the peak is relatively small with respect to that 
of K-' and V'; this is just the hallmark of a competition 
between the two contributions. 

At a higher temperature T ~ located in between the 
two peaks of Cy, IC' has a flat minimum and V' a flat 
maximum. Finally, at still higher values of temperature, 
]C' exhibits a second maximum at T3, and V' is smoothly 
decreasing; in correspondence, Cy exhibits the second 
peak, of kinetic origin. 

As U is increased (fig.(b)), the value of Ti decreases and 
the absolute height of both the above contributions dras- 
tically vanishes, so that the low-temperature peak be- 
comes a sort of 'shoulder'. At the same time, the mini- 
mum of the kinetic contribution and the maximum of the 
potential contribution located around T2 have become 
more pronounced, and T2 itself has decreased (see fig.(c)). 
As U reaches the value 2t of the metal-insulator transi- 
tion, both Ti and T2 vanish and the magnitude of the cor- 
responding extrema become infinite (fig.(d)). For U > 2t 
(figs, (d) and (e)), the T2 extrema are regularly restored 
and, since the Ti-extrema have disappeared, they become 
the new low-temperature extrema. At this temperature 
Cv exhibits now a new peak. Thus, for U > 2t, the po- 
tential contribution prevails on the kinetic one, and the 
nature of the low-temperature peak changes with respect 
to the case U < 2t. Notice that T2 now increases with U 
(figs.(e) and (f)). Finally, at higher temperatures another 
broad peak originates from the (old) second maximum of 
the kinetic part. Such high- temperature (T3) peak is very 
broad, and it definitely disappears when U is further in- 
creased above 3i. 

The above observations show that at half-filling, passing 
through the point U = 2t oi the metal-insulator tran- 
sition, the nature of the low-temperature peak changes 
its origin from a kinetic to a potential one, whereas at 
moderate J7's a further peak of kinetic origin appears 
at higher temperatures. In pass let us also notice that 



at strong coupling a kinetic (potential) peak is a peak to 
which only IC'iV') basically contributes, V'{1C') being al- 
most vanishing (see fig.^^; in contrast, at moderate J7's 
a kinetic (potential) peak is a peak for which the kinetic 
contribution slightly prevailing on the potential(kinetic) 



one. 
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FIG. 15. The temperature dependence of the kinetic 
(dashed) and potential (dotted) contributions to the spe- 
cific heat (solid), at half filling and different moderate val- 
ues of U . Contrary to the case f/ 3> t of figjl4| at mod- 
erate coupfing /C' and V' are competing, since they have 
relatively large contributions of opposite signs at roughly 
the same temperature. The peaks of Cv are thus 'kinetic' 
(resp. 'potential') when K,' (resp.P') prevails on the other. No- 
tice that the low-temperature peak changes its origin from ki- 
netic to potential across the metal insulator transition point; 
the high-temperature one, present up to U /t ^ 2.5, is instead 
always of kinetic origin. (For editing reasons the two bottom 
figures have a different y-axis scale.) 

The results obtained for our model can be compared 
with those concerning the ordinary Hubbard model. 
In the strong coupling regime of this model the low- 
temperature peak is attributed to spin excitations (the 
corresponding temperature being of the order of J = 
4t^/J7), whereas the high-temperature peak is related to 
the charge excitations (since it is located at fc^T ~ U). 
With lowering U, it is widely accepted that the two peaks 
merge at f/ ~ 4t; however, some investigations have been 
carried out at still lower J7's, showing that a double peak 
structure reappears for D=l [|ll| and T)—2 [|l6|,0. It is 
customary to relate the origin of these new peaks again 



12 



to spin and charge degrees of freedom respectively. 
The Hubbard model is considered the paradigm within 
strongly correlated systems, so that the presence of a 
two-peak structure in the specific heat of such systems 
tends naturally to be interpreted as the signature of spin 
and charge excitations. 

However, in the authors' opinion, not enough attention 
has been devoted to the effect that further interaction 
terms in the Hamiltonian have on the specific heat. To 
this purpose, the exact results obtained for our model 
show that, when a possible competition between single 
and paired carriers is taken into account, the specific heat 
turns out to exhibit a structured two-peak behavior, in 
spite of the fact that only charge degrees of freedom are 
involved. Although our model neglects some terms such 
as nearest neighbors charge interaction Vhiafijui)^ we 
believe that it can reproduce some features of realistic 
materials which are not explicitly taken into account in 
the ordinary Hubbard model, namely: a) the opening of 
the gap at a finite value of U/t^ i.e. at a finite value 
of pressure on the sample; b) the lack of particle-hole 
symmetry, observed in heavy fermion compounds; c) the 
presence of a mechanism favoring the kinetic of paired 
carriers, as it is the case in cuprate superconductors. 
In view of these observations, we suggest that the inter- 
pretation of a two-peak structure in Cv may not neces- 
sarily be related to spin and charge excitations; a com- 
parison with the behavior of pure spin quantities, such 
as magnetic susceptibility, in correspondence of the peaks 
temperature would be more probative. 



VI. CONCLUSIONS 



In this paper we have calculated the exact thermo- 
dynamics of an Extended Hubbard model by means of 
Sutherland Species technique, which we had previously 
introduced to determine the ground state properties of 
the same model [|]. The model describes a competi- 
tion between the dynamics of single carriers and that of 
short-radius paired carriers; such competition is modu- 
lated by the values of the electron filling p and the on-site 
Coulomb repulsion U . We have calculated the partition 
function of the model, and derived the finite tempera- 
ture behavior of different physical quantities, namely the 
chemical potential, the compressibility, the local mag- 
netic moment, and the specific heat. We have discussed 
the changes of such observables across the point of the 
metal-insulator transition U = 2t, providing explicit low- 
temperature expressions for Cv and /i; in particular fi is 
found to undergo an unusual transition from a linear to a 
TlnT dependence. We have then focused on the specific 
heat, which turns out to exhibit interesting features, such 
as a nearly universal crossing point and a double peak 
structure. The two peaks, which are shown to be related 
to charge degrees of freedom only, are present in ranges 
of U/t both below and above the metal-insulator tran- 



sition value. We have discussed the two peaks in terms 
of kinetic and potential contribution to the spectrum, 
outlining the differences between the cases of strong cou- 
pling and moderate coupling, and comparing our results 
with that of the ordinary Hubbard model. 



The method presented here to derive the partition 
function of our model can be applied, with straightfor- 
ward generalization, to further integrable extended Hub- 
bard models j26j involving two Sutherland species. Work 
is in progress along these lines. 
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